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Abstract 

We review the non linear sigma model approach (NLSM) to spin chains and spin lad- 
ders, presenting new results. The generalization of the Haldane's map to ladders in the 
Hamiltonian approach, give rise to different values of the parameter depending on the 
spin S, the number of legs and the choice of blocks needed to built up the NLSM fields. 
For rectangular blocks we obtain = or IttS depending on wether n^, is even or odd, 
while for diagonal blocks we obtain = lirSrif. Both results agree modulo 2tt, and yield 
the same prediction, namely that even ( resp. odd) ladders are gapped (resp. gapless). 
For even leeged ladders we show that the spin gap collapses exponentially with ni and we 
propose a finite size correction to the gap formula recently derived by Chakravarty using 
the 2+1 NSLM, which gives a good fit of numerical results. We show the existence of a 
Haldane phase in the two legged ladder using diagonal blocks and finally we consider the 
phase diagram of dimerized ladders. 

* Based on a talk delivered at the Summer School on "Strongly Correlated Magnetic and Superconductin 
Systems", held in Madrid, Spain, July 1996 







Introduction 



There have been three major developments in the 80's and 90's in Condensed Matter or more 
specifically in Strongly Correlated Systems. In historical order they are 

• Haldane's conjecture in Id antiferromagnetic spin chains (1983) [|l], H 

• Discovery of high-T^ superconductivity and antiferromagnetism in doped and undoped 
cuprate compounds (1986) [H 

• Discovery of Ladder Materials (1987,91) 0, § 

All these findings have in common various features: low dimensionality, antiferromagnetism, 
importance of quantum fluctuations and the fact that they constituted theoretical and exper- 
imental surprises. Moreover these topics can be studied using sigma model techniques, which 
establish a methodological link between them. In this talk we shall be mainly concerned with 
spin chains and ladders. 

The first surprise deals with the behaviour of Id antiferromagnetic Heisenberg spin chains 
(AFH) as a function of the spin. The Hamiltonian describing the AFH model is given by, 

L 

E,^^^ = jY.'^{n)^{n^\) (1) 

n=l 

where J > is the exchange coupling constant and S(n) is a spin-S matrix acting at the n*'' 
site of the chain. Haldane's conjecture concerning the spectrum and spin correlations of the 
Hamiltonian (|I]) is given in table 1 , 



2S 


Spectrum 


Correlations 


even 


gapped 


exponential decay 


odd 


gapless 


algebraic decay 



Table 1: Behaviour of spin chains 



Nowadays there is sufficient theoretical, numerical and experimental evidence to support 
this conjecture, which should rather be called a "theorem" despite of the lack of a rigorous 
mathematical proof 0. 

The third surprise, in the historical order shown above, deals with the behaviour of spin 
ladders as function of the number of legs ( for a review and references on the subject see 

1^ ) . A ladder is an array of spin chains coupled as in figure 1 . The Hamiltonian of the system 
is given by. 



H,e, = JYTaU ELi Sa(n)S,(n + 1) (2) 

where Sa(n) are spin-S matrices located in the a^^ leg at the position n = 1, . . . , L. We consider 
periodic boundary conditions alongs the legs, i.e. Sa(n) = Sa(n + L), and open EC's along the 
rungs. The intraleg coupling constant J is positive but the interleg coupling constant J' can be 
either positive or negative. The qualitative behaviour of spin 1/2 ladders is given in table 2, 
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Spectrum 


Correlations 


even 


gapped 


exponential decay 


odd 


gapless 


algebraic decay 



Table 2: Behaviour of Spin 1/2 Ladders 

This result holds for both signs of J'. The analogy between the integer/half- integer behaviour 
of spin chains and the even/odd behaviour of spin 1/2 ladders is evident. In the limit where 
J' goes to minus infinity, the spin ladder Hamiltonian (Q) becomes equivalent to the spin chain 
Hamiltonian (|ip for an effective spin equal to rig /2, and hence the ladder behaviour given in table 
2 follows from the behaviour of the spin chains. What is not so obvious is that this behaviour 
holds not just for strong ferromagnetic rung couplings but also for antiferromagnetic ones, re- 
gardless their magnitude. As for the Haldane conjecture there is by now sufficient evidence 
to support the "ladder conjecture" coming from Quantum Monte Carlo, exact diagonahzation, 
mean field theory, experiments on ladder materials like VOPO or cuprates, bosonization, finite 
size analysis, sigma model, etc pf. 

In this talk we shall give a unified description of the behaviour of spin chains and spin ladders 
utilizing the non linear sigma model (NLSM) |10], 0. Let us first recall the logic underlying 
Haldane conjecture. It is based in the following map|I|, 0, 



Low Energy Modes of the Spin Chain — > Non Linear Sigma Model (3) 
which is obtained in the semiclassical limit where the spin S becomes very large, although @ 



can be derived on more general grounds based on symmetry arguments [|T^. However these type 
of arguments miss the theta term in the action, which plays a crucial role in determining the 
physics of the model. 

Hence using the properties of the (NLSM) one derives those of the low energy spectrum of 
the spin chain. 

To arrive at the results of table 2 for the ladders we shall follow the same logic as for spin 
chains, that is, we shall construct a map 



Low Energy Modes of the Spin Ladder 



Non Linear Sigma Model 



(4) 



and use the properties of the NLSM. In this way the NLSM provides a unified and economic 
approach to both problems. Actually, the 2d Antiferromagnetic Heisenberg model, can also be 
mapped into the 0(3) NLSM in 2+1 dimensions |l^. The RG analysis at finite temperature 
of the later model made in has unravelled many of the properties of the undoped cuprate 
compounds ( second surprise). On the other hand the NLSM's that appear in the r.h.s. of 
the maps (H) and (H) are defined in 1+1 dimensions. This poses the problem of the crossover 



between Id and 2d Heisenberg systems through spin ladders [15 . 



The Non Linear Sigma Model: A Primer 

Before we construct the maps (^ and (^, it will be convenient to review the basics of the 0(3) 
NLSM. This model was proposed as a toy version of QCD, as can be seen by enumerating its main 
features: asymptotic freedom, dynamical mass generation, instantons, skyrmions, integrability , 



2 



etc (for a review see |16, [T^])- For these and other reasons the NLSM still attract the interest of 
many physicists and mathematicians. In the case of spin chains and ladders the NLSM becomes 
not just a toy model but a reallistic model describing the low energy physics!. 

The 0(3) NLSM is a relativistic quantum field theory in 1+1 dimensions whose field $ is a 
three component vector living on the 2d- sphere S'^, 



$2 = 1, 

The euclidean action of the model is given by, 



S 



d X 



{d,^f + i^ e^, $ ■ (a^$ X d,^) 
Zg oTT 



(5) 



(6) 



where 5^ > is the sigma model coupling constant, and e^;^ the 2d Levi-Civita symbol. The 
quantity 



W = — d'xe,,^ 



(9^$ X d,^) e z 



(7) 



takes an integer value for field configurations $(x) which go to a fixed value, say $0, when 
|x| ^ 00 ( this condition is required to have a finite action). Compactifying the space-time 
into the sphere S'^, the integral (0) gives the winding number of the map S'^ (space-time) — > 5*^ 
(target space). The parameter 9 enters in the partition function as e'^^, and therefore the 
integer character of W implies that 6 is defined modulo 27r. The theta term in the action is a 
truly topological term, which leads to dramatic non perturbative effects. 

Notice that W changes its sign under a parity ( or time reversal) transformations. Hence 
the topological term of the action breaks explicitely parity (or time reversal) unless = or vr. 
Indeed if ^ = (mod 2 tt) the topological term is completely absent while if ^ = vr the winding 
number contributes to the action with a sign, i.e. (—1)^. We expect from these properties that 
the AFH-spin chains and ladders which are parity invariant will be associated with ^ = or vr. 

The basic properties of the NLSM for these values of 9 are given in table 3 |18|, |l^, [|, 



6 ( mod 27r) 


Spectrum 


Correlations 





gapped 


exponential decay 


TT 


gapless 


algebraic decay 



Table 3: Behaviour of the 0(3) NLSM 

This behaviour, which is independent of the magnitude of the coupling constant g, will allow 
us to make a quick derivation of table 3 in the strong coupling limit. Before we do that, let's 
write the Hamiltonian that follows from the action (BI), 



NLSM 



dx 



9 ^ 



e 

4:71 
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where $'(x) = dx^{x), c is the "speed of light" and \{x) is the angular momentum density, 
which is defined as, 



1 = $ X — — 
dt 



(9) 
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Besides the constraint the fields $ and 1 satisfy, 



l(x)-$(x) = (10) 
and the cannonical equal-time commutation relations, 

[/'^(x), (''{y)] = ie^Hix - y)l''{x) 

[/'^(x), ^\y)] = ie"^^5(x - (11) 
[<l>"(x), $''(?/)] = 

It is very illustrative to prove the statements given in table 3 in the strong coupling limit 



((7 >> 1 ), for which we shall use a regularized lattice version of the Hamiltonian (§) pO|, ^ 



H'nLsIi = ^ E U 1' - 7 ^n^n+l ) (12) 




At every site of the lattice there is a "particle" moving on a sphere ($^ = 1), with angular 
momenta In- If ^ = the angular momenta takes all possible integer values. In = 0, 1, . . . cxd 
pi| . However if 6 = it there is a monopole of charge 1 at the center of every sphere, which 
implies that the possible values of the angular momentum are restricted to half-integers values, 
/„ = l/2,3/2,...oo m. 

In the limit g » 1 the kinetic term g J2n dominates over the potential term ~^ ^n^n+i, 
and the ground state is obtained choosing the smallest possible value of /„ at every site. 

If ^ = the optimal choice is given by In = 0, Vn, which yields a unique ground state. The 
first excited states are obtained by choosing the irrep / = 1 in one site and Z = in the rest 
of the chain. Since this can be done at every site, there is a huge degeneracy, which is broken 
by the potential term, which delocalizes the / = 1 excitations. The 3L degenerate first excited 
states become a band of Z = 1 magnons, separated from the ground state by a gap, which can 
be computed in perturbation theory. In this gap was computed up to 6^^ order in 1/g'^. The 
first three terms read, 

A = cg(l-~ + O.O74I + 0(1)) (13) 

In the weak coupling limit {g << 1) a. perturbative RG analysis shows that the gap vanishes 
exponentially as [|TB| 

A ~ -e-^-/9 (14) 
9 



The proportionality constant depends on the regularization of the model 12^ (see also 



Combining the strong and weak coupling analysis and using Fade aproximants the authors of 
pO| concluded that the gap should never vanishes as long as g is non zero. 

For 6 = TT the kinetic energy is minimized by the choice /„ = 1/2 Vn, which still leaves a huge 
degeneracy for the ground state. This degeneracy is lifted by the potential term, which leads to 
an effective AFH model when restricted to the subspace In = 1/2 Since the S' = 1/2 AFH 
model is massless one gets that the NLSM at ^ = vr is also massless. In reference it was 
argued that this gapless behaviour persists to all values of g. 
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Haldane map for spin chains 



The map from the Heisenberg model into the NLSM can be done in various ways: using coherent 
states in the path integral formalism, generalizing the Hubbard-Stratonovich formula in the 
partition function, or applying gradient expansions in the Hamiltonian formalism. We shall 
mainly use the later one 0, but we shall also briefly explain the path integral approach for 
ladders. 

The starting point of the construction is the spin wave analysis of the AFH model. This 
consist in the linearization of the evolution equations of the spin operators S(n), around the 
classical Neel configuration. The equations of motion of the spin operators read. 



^ = i[i/,hain,SH] (15) 

= - JS(n) X [S{n + 1) + S(n - 1)] 

The basic asumption is that S{n) deviates by a small amount from the alternating Neel 
configuration, 

S(r2) = (-l)"5z + s(n) (16) 

where S is the magnitude of the spin and z is the unit vector in the vertical direction. In the 
linearized approximation eq.(|l5]) becomes, 

^ = i(-l)"+V^[C(n + 1) + C(^ - 1) + 2C(n)] (17) 
where C{n) = s^{n) + is^(n). The spin waves are the plane wave solutions of eq.(|T 



({n) = e'('^*+*^") (tPik) + (-1)"+ V(^)) (1^ 
We are interested in the low energy and long wavelength solutions which are given by. 



uj = vk 

ip{k) ~ Ak (19) 
(f){k) ~ B 

Equation (|T^ is satisfied by the ansatz (0) if 2A = B. The spin wave velocity v is given by. 



V = 2JS (20) 

The two transverse spin wave solutions (p^ can be identified with the massless goldstone 
bosons associated with the breaking of the rotational symmetry of the AFH Hamiltonian by the 
classical Neel state. Eq(|TB|) implies that the spin variables have two slowly varying components 
around the momenta k = { field 1) and k = n { field $), which correspond to the local spin 
density and the local staggered magnetization respectively. 

This elementary spin wave analysis is very useful in order to construct the map from the 
Heisenberg model into the NLSM. This map can be mathematically formulated as a change of 
variables 0. First of all, one has to divide the chain into blocks of two sites ( see figure 2). The 
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block, say (2n, 2n+l), is given the coordinate x= 2n+l/2 and for every block one performs the 
following change of variables, 



S(2n) = l(x) - S^{x) (21) 
S{2n + 1) = \{x) + S^{x) 



The inverse of eqs.(pl|) are 



= |[S(2n + l) + S(2n)] (22) 
$(x) = ^[S(2n + l)-S(2n)] 

These relations imply that 1 is a local spin density and $ is a local staggered magnetization. 
Using ( [21| ) and the eqs S^(n) = S{S + 1) one gets, 

\\x) + S^^^{x)=S{S + l) (23) 
1- $ = 



In the semiclassical limit S >> 1 eqs.(^) become the sigma model conditions (^ and (pl)|). 
The commutators (|TT|) can also be obtained in the limit S >> 1 from the commutation relations 
of the spin matrices, with the following identification of the Dirac's delta function and the 
Kronecker's delta symbol in the limit where the lattice spacing 6 goes to zero, 

Six-y) = \im^ 6,,y (24) 

The term 26 in the denominator of eq. (|2^) is the lattice spacing in the x- variables ( from 
now on we shall set 6 = 1). 

The AFH Hamiltonian (||) reads in terms of the new variables 1 and 

i^chain = E.[-S{S + 1) + 2\^{x) + \{x) \{x + 2) (25) 
+S{^{x) \{x + 2) - l(x) $(x + 2)) - S^^{x) $(x + 2)] 

Making the asumption that the fields l(a;) and $(x) vary slowly in the x coordinate, we can 
perform a gradient expansion and truncation of the higher derivative terms. Keeping terms with 
at most two spatial or time derivatives in the field we obtain the following Hamiltonian, 

^(Continuum) ^J^^ ^^^2 _ ^(j^/ ^ ^/j) ^ S^^'f] (26) 

which coincides with the NLSM Hamiltonian given in (^) upon the identifications, 

e = 27fS, g = ^, c = 2JS (27) 

This is the desired result which permits to derive table 1 from table 3. From eqs (0) and 
(P7|) we can estimate the value of the spin gap and the correlation length as functions of the 
spin. 
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ttS 



S 



(2^ 



In table 4 we show the numerical values of the gap, correlation length and spin velocity of 
the S= 1 and 2 Heisenberg chains, obtained using Quantum Monte Carlo |24, ^ and DMRG 
methods |2^. 



Spin 


A/J 




c/J 


S = l 


0.4107 


~ 6 


2.47 


S = 2 


0.049-0.085 


49 


4.16 



Table 4: Spin gap of the S' = 1 and 2 chains. 

This table shows the semi-quantitative agreement between the numerical (exact) results and 
the NLSM estimates In the semiclassical limit where S —>■ oo the NLSM predicts that the 
gap should go to zero exponentially fast, recovering in that way the classical gapless behaviour. 
This shows that the existence of the gap is a truly quantum mechanical effect. 

To end up our review of the spin chains, we would like to make some comments on the non 
uniqueness of the map (|]). Indeed eqs(^) and which give the map between the spin and 
the sigma model variables, depend on the choice of the 2 site blocks. The other possible choice, 
given by the blocks (2n+l, 2n+2), leads two another couple of variables 1 and which are 
linearly related to 1 and by the eqs. 



1 = 1-5$' (29) 
$ = $ 

These eqs. leave invariant the constraints (^), ([T0|) , and in fact they can be obtained by the 
following transformation 0, 

J ^—iS J dx fj' cos« J ^iS ^ dx 13' cosa (30) 

^ ^—iS J dx P' cosa ^ ^iS J dx f3' cosa 

where a and /? are the spherical coordinates that parametrize the staggered field 

$ = (sina cos/3, sina sin/3, cosa) (31) 



However the Hamiltonian (p6D is not left invariant, but the change affects only the theta 



parameter which for the new variables becomes, 

e = -2718 (32) 

This change in the value of 6', depending on the blocking, can also be seen as due to a 
parity transformation. Recall that parity changes the sign of the topological term W. The two 
blockings (2n, 2n +1) and (2n+l, 2n+2) are indeed related by parity. In the next section we 
shall consider more general blockings which will lead to changes in 9 but that are not related 
by parity transformations. If 6' = vr the change ( P^ ) leaves invariant the sign factor (— 1)^ 
appearing in the partition function. 
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Haldane map for spin ladders 



There are two ways to study the problem of spin ladders using NLSM methods. The first one 
consists in the application of the Haldane's map to every chain forming the ladder, obtaining 
a system of sigma model fields coupled by rung interactions. This approach is similar in 
spirit to the bosonization studies of spin ladders made in references , and its applicability is 
reasonable in the weak coupling regime J' / J << 1. The other possible approach is to built up a 
unique sigma model field describing the low energy modes of the spin ladder as a whole [1^, |TT|. 
This approach should be appropiate to study the intermediate coupling regime {J' / J ~ 0(1)). 
The strong coupling regime (J'/ J >> 1) has been mainly studied using perturbative [|^, ^ ISTJ 
and mean field methods ||32|. We shall see that the later approaches are related to the discrete 
version of the NLSM based on the Hamiltonian ([T2|) . 

The map from the ladder into the NLSM follows the same steps as that of the spin chain. 
We shall review below the approach of reference [lin . 



Spin Wave Analysis of Ladders 

This analysis shows that at the linearized level there are two massless modes corresponding 
to two Goldstone bosons. These modes describe the gapless deviations from the classical Neel 
solution, which in the case of AF-coupling along the rungs reads, 

S'^'^'in) = (-l)°+"5 (33) 

In the rest of the talk we shall confine ourselves to this case. The transverse components of 
the spin waves , i.e. Ca{n) = s^{n) + is^(n) , are given by 

(^(n) = e'('^*+'=") {Aak + {-ly+^'+^B) (34) 

where u = vk. The quantities Aa represent the fraction of total spin carried out by the a*''— leg, 
and we shall normalize them as, 

E^a = l (35) 

a 

The spin wave velocity v and the A'^s are given by. 



{^y = Jne/j:a,bL-l (36) 

= Y.b La},/ Y.c,d -^ci 

where is the inverse of a matrix L defined as follows, 

{4 J + J' a = b = loTHe 
AJ + 2J' l<a = b<ng (37) 
J' |a - 6| = 1 

Besides the 2 Goldstone bosons there are, at the linearized level, 2{ne — 1) massive modes. 
As an example we give below the values of their masses for the ni = 2 and 3 ladders. 
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m„,=2 = SVSJT (3^ 
m„,=3 = SJj'{J + AJ') , 5 Jj'(J+ 12 J') 



Of course in the hmit when J' — > the masses of the massive modes (|38| ) dissapear and all 
the modes become massless, as should be the case for a set of uncoupled chains. The basic 
asumption in the construction of |TD| is the truncation of the massive modes, keeping only the 



gap less ones, which is justified if the energy scales, generated non perturbatively, are lower than 
the gap of the massive modes. One may expect that this condition is satisfied in the intermediate 
and strong coupling regimes. 



Map: Spin Ladders ^ NLSM 

The spin wave analysis serves as a preparation for the more complicated job of finding the map 
from the spin ladder into the sigma model. In the Hamiltonian formulation the first step is to 
split the ladder into blocks of sites, defining for each block an average angular momentum 1 
and staggered magnetization $ as follows. 



l(x) = ^E(a,n)eB(.) S,(n) (39) 

where x denotes the center of mass position of the block B{x) along the leg axis, and the 
prefactor ni/uB is required in order that 1 satisfy the cannonical commutation relations ([TT|) in 
the continuum limit. The relation (p^ ) reads in the more general case. 

For single chains there are only two types of blockings related by parity. However for ladders 
there are many different types of blockings, which may in principle lead to different results. 
In order to choose those blockings that are physically acceptable we shall impose the following 
conditions 

• The blocks must have an even number of sites. 

• Every block must contain a site (or more ) of every leg the same number of times, which 
implies that should be a multiple of n^. 

• Every site in a block must have a nearest neighbour belonging to the same block. 

Motivated by the spin wave solution (0), we shall propose the following ansatz for the spin 
operators in terms of the sigma model fields, 

Sa{n) = Aa \{x) + (-1)"+"^ $(x) for {a,n) e B{x) (41) 
The consistency between eqs. (|39|) and (^T]) is guaranteed by the following identities, 
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I](a,n)eB(a;)( — 1)"^" — 

E(a,n)eB{aO = (42) 

I](a,n)eB{a;)(~l)" "^a = 

which can be proved from the conditions on the allowed blocks given above and the symmetry 
properties of Aa. 

Strictely speaking we should add to the r.h.s. of ( ^T]) a set of fields describing the massive 
modes whose existence we discussed in the spin wave analysis. Their inclusion makes the map 
(PI) more rigorous from a mathematical point of view, allowing a careful derivation of the 



effective Hamiltonian governing the dynamics of the fields 1 and $ [1^]. Having found the 



effective Hamiltonian the massive fields are discarded, so that we shall not consider them from 
now on. 

Let us now consider a few examples ladder's blockings. 

Columnar Blocks {ns = ng) For ladders with an even numbers of legs the smallest possible 
blocks, satisfying the conditions given above, coincide with the rungs (see Fig 3), i.e. ub = rii. 
The partition of the ladders into rungs is what one effectively does in the study of the strong 



coupling limit of ladders ^ . In that sense one may expect to find relations between the 

NSLM approach and the strong coupling analysis. 

Replacing the ansatz (^Tj) into the ladder Hamiltonian (^), and taking the continuum limit 
we get the following expression. 



Comparing (|^) with the NLSM Hamiltonian (^ we find. 



(43) 



e= 

s 



'^'^ Jn./EafiL-l (44) 



9= 1/ (S^JniY.a,b La 



Since 6' = we see that the NLSM is gapped, which implies that even ladders should always 
be gapped for any value of the spin. The velocity c coincides with the spin wave velocity v (^6]). 
Finally, the coupling constant g has a non trivial dependence of the number of legs and the ratio 



J' / J. This is interesting because we can derive from the eq.(0) the dependence of the spin gap 
on the ladder's parameters. 

For odd ladders the minimal blocks satisfying the conditions above must have at least ub = 
2ni sites. Unlike the blocks with ub = ne there are now more geometries. We shall in what 
follows consider the cases of even and odd values of n£. In figures 4 and 5 we show two of them. 
The choice of figure 4 was the one studied in [|TU|. We shall give for completeness the results 



obtained there. Later on we shall study the blocks of figure 5. 

Rectangular Blocks {ub = 2nf) The Hamiltonian that one obtains after taking the 
continuum limit is (fig 4), 
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H 



dx 

Y 



a,b a=l 

In it was shown that the value of 6 corresponding to (^5]) is given by, 

e = 



(45) 



n£ : even 
2718 rii : odd 



(46) 



The vanishing of 6 for even ladders follows simply from the symmetry property Aa = A^^+i-a, 
while the value 6 = 27:3 for odd ladders requires some more work |Ty]. The other NLSM 
parameters of (|4^) are given by. 



Jne 



Z^a.b a.b 



fay ,L-i 



9— V {S^2JniJ2a,bLal - ISniMd) 



(47) 



where 5n^,odd = 1 ( resp. 0) if ne is odd ( resp. even). For ne even the delta terms appearing 
in these eqs. are absent and we recover (^^, except for a renormalization of c and g (c^ect = 

Ccol, greet = fl'col/\/2)- 

For rii odd, the formula ( ^Tj) does not coincide with the spin wave velocity (p6|), except for 
the case of the spin chains ( = 1)!. For large ladders the velocity approaches the value 475* ( 
J' = J) which differs by a factor of ^/2 with respect to the 2d value which is 2\/2JS. 



Diagonal Blocks [ub = 2ni) Within the gradient approximation which we are using, it 
is consistent to assign to all the points within a diagonal block (fig. 5) the same coordinate x. 
Under this asumption we get the following Hamiltonian, 




where pa is defined as 

r J + J'/2 a = lorn, 
~ \ J + J' else 

The value of 6 that comes out from ( ^5]) is given by, 

= 87r^ ^"/°^° = 87.3 VaL-l (50) 

l^b,c ^b,c ^b^c a,b 

To evaluate 6 it is convenient to write (^Up in matrix form as follows, 

9 = 87TSne < F\PL-'^\F > (51) 

where \F > denotes a normalized n^-component vector with all its entries equal to 1/y/ne and 
P is a diagonal matrix with its entries given by pa- The trick is now to write the matrix L in 
the form L = 4P — K~ , where K~ is defined as. 
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1 a = b = 1 or n£ 
K-^ = j'y^^ 2 a = b ^1 or rie (52) 
-1 |a-6| = l 



Making a "Dyson" decomposition of L ^ 



4P 4P^~4P-K- ^^^^ 



and using the fact that K annihilates the vector |F > we finally get, 



e = 2nSni (54) 



It is quite remarkable that all the coupling constants in the expression ( pOD just combined in 
order to produce a result which only depends on "global" data S and ng. Eq.(|51D agrees with 
(^) mod 27T, implying that the choice of block does not affect this parameter. Let us see what 
are the values of g and c for the Hamiltonian ( |48|), 



(i) 



2 



2 

^ 1,6 \ 



(55) 



Path Integral Derivation of Haldane Map for Ladders 

In this section we shall apply coherent state techniques to derive the map from the spin ladders 
into the NLSM (for a review see reference [^). The map for the 2 legged ladder was first worked 
out in 0. Our results agree with the result of this reference and extend them to generic values 



of ( see also |TT| for an alternative derivation of the results given below). 

The action of the spin ladder system in the path integral formulation can be written in the 
following form, 



5 = / rfr {i E„ TTaU AiSain)) ■ (56) 

-JEnEaLl Sa{n)Sa{n + 1) - J' EnEa=l^ (n) S^+i (ri) } 

where Sa(n) is a classical spin variable satisfying Sa(^)^ = 5'^ and A(Sa(n)) is the vector 
potential that fuUfills the constraint rotA(n) = n (n = S/S*). 

The long wavelength limit of the ladder's spin variables Sa(n) is given, according to (0), by 

S,(n) = SAJin) + {-ir+-S^{n) (^1 - ^^A'f^ ' (57) 

where the term with the square root is needed for the correct normalization of the variable 
Sa(n). For spin chains this later term can be replaced by 1, but for ladders it gives a non trivial 
contribution when considering the coupling between rungs. 

Introducing ([57| ) into (^) and performing the standard gradient expansions one finds the 
following NLSM parameters 
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6 — 2txS 5ni,odA 



2 

5 



9= 1/ [S^JniJ2a,b La 



which coincide with those obtained using columnar blocks (|4^) for even legged ladders. For odd 
ladders the results of the path integral and those obtained with rectangular blocks differ, except 
in the case of spin chains {ne = 1). This poses the problem between the relation between the 
Hamiltonian formalism and the path integral formalism for this type of ladders. In any case, 
notice that the path integral approach gives for even and odd legged ladders a value of c identical 
to the spin wave velocity (p^). 

Before we extract more consequences from the mapping of the ladders into the NLSM, it is 
convenient to express c and g in terms of 



The function 

Let us define as |10| 



rii 



J2a,b La 



(59) 



n^jOdd 



m=l,3,...,ni—l 



sm 



fe) + * 



cos 



2 nm 
2ni 



Its explicit expression and numerical values for z = 1 m the cases = 1, 2 ,3, 4 and oo is given 
in table 5. 





1 


2 


3 


4 


oo 




1 


1/(1 + 1) 


(l + ^)/(l + f) 


(l + f)/(l + ^ + f) 


1 

1 + 2 


/n.(l) 


1 


0.666 


0.6190 


0.5882 


0.5 


2 + „,v^) 


0.853 


0.677 


0.6178 


0.5884 


0.5 



Table 5 

fnt{z) is a monotonically decreasing function in which varies from 1 to 1/n^ 5nf,odd as z 
varies from to oo. Interesting enough, the sum appearing in eq. (p9|) can be performed yielding. 



+ Z 



1 + f (i + VTT 



ni 



n, (1 + zf^ (l + ! (l + Vl + t)Y' + (-1)"^ 



In the isotropic case {z = 1) we get. 



(60) 



■'"'^ ' 2 i^' + „,V2 (3 + 2V2) ' + (-!)"<, 
For > 2, a good approximation of (^) is given by ( see table 5), 



(61) 
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1 
2 



1 



V2, 



(62) 



In eqs ( |60|) and (^) the difference between even and odd ladders dissapears exponentially 
as rii increases. 

This ends the review of the properties of and we return now to our general discussion. 



Summary of Results and Conclusions 

In table 6 we summarize the results obtained so far for the three types of blocks we have 
considered so far, 



Block 


ni 


9 


c 


9 


Columnar 


even 







2/Sni,^fn, 


Rectangular 


even 







^/2/Sni^fnt 


Rectangular 


odd 


27tS 


/ \ 1/2 
2V2JS / 1 1 \ 


/ \ -1/2 
^ (l 






Diagonal 


both 


27TSni 


g[2/„,(l + (l-i,)^)-l]^^^ 




Path Int. 


both 


2nS6ni,odd 


2JS/^h^ 


2/Sni,^fn, 



Table 6 



From this table we can extract the following conclusions, 

• The values of 6*, c and g are block dependent. For example, for even we get g^^ct = 5'coi- 
Since the size of both blocks is different we could interpreted the previous relation as the 
RG relation bewtween g's at two different length scales (ns(rect)/nB(col) = 2). Curiously 
enough the value of g^ec is smaller than the value of 5^001, suggesting that this "one step RG 
flow" is similar to the RG flow of the 2+1 NLSM, where the coupling constant g decreases 
at longer distances in the " renormalized classical region" [Q. After this first RG step, 
which truncates the ladder's degrees of freedom down to those of the 1+1 NLSM, the 
value of g will start growing as a consequence of the 1+1 NLSM RG equations, so that the 
physics of the system will be dominated by the strong coupling regime. 

From table 6 we find two " RG-invariant" quantities: the 9 parameter and the perpendicular 
spin susceptibility. 

• The invariance of the 9 parameter relies on its periodicity, which implies that for all 
purposes we may take 9 = 27iSni. This is really a topological result for it only depends 
on the "global data" S and rii, and not on the values of the ladder coupling constants J 
and J'. For odd ladders this result must be a consequence of the generalization, due to 



Affleck [^, of the well known theorem by Lieb-Schultz-Mattis [0, that asserts that in 
the infinite length limit the odd ladders must either have a degenerate ground state or else 
there are gapless excitations. 
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• The value of the bare perpendicular susceptibility x° is given for all the block choices by, 

Notice that for large nn the susceptibihty per site jni goes to a finite value. 
It is quite interesting to compare (O) and the spin wave velocity 



2.JS , , 

(64) 



with the corresponding expressions of the bare perpendicular spin susceptibility and spin wave 
velocities of a dimensional sigma model | 37| , 

where a denotes the lattice spacing. 

The comparison of (0) , (|6^ and (|65|) suggest somehow that ladders behave as (iiadder dimen- 
sional spin systems with, 

(^ladder = (66) 

This naive definition of "fractal" dimension of ladders helps to explain some numerical facts. 
First of all, if we choose J' = 0, then /n^(0) = 1 and we get c^iadder = 1, which indeed corresponds 
to Id chains. For the isotropic models and large we get, 

2 

(^ladder ~ ^ i (67) 

which converges towards d = 2 (plane) from below. 

A less heuristic proposal is to associate with a "finite size" effect in the renormalization 
constants Zc and of the ladder. Combining the recent work of Chakravarty [|15|, together with 
the results presented above, we shall propose a finite size correction of the spin susceptibility, 
spin velocity and spin-stiffness renormalization constants of ladders as follows, 

Z^{S,ne) = 2fn,Z^{S) 

Z,iS,ne)= Ze(S)/y2/~ (68) 
ZpAS,ni)= Zp^{S) 

where Z^{S), Zc{S) and Zp^S) are those of the 2d spin system (i.e. tit = oo). The last eq. in 
(|68|) follows from the first two thanks to the relation Zp^ = Z^Z^. We shall give some numerical 
support to (p8| ) in the next section. 

The general formulation we have introduced above, will allow us to discuss certain important 
issues concerning the ladders. 
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Spin Gap of Even Ladders 



Choosing the columnar-block description of the even ladders we deduce from table 6 and eq. ([14D 
the following expression for the spin gap, 



Al^+^) ~ JS"^ Uf, exp i-T^Sn 



(69) 



which predicts a exponential decay of the g function on the number of legs [|T^]. This 

implies in particular that the spin gaps of the 2, 4 and 6 legged ladders should be related. Indeed 
one finds from numerical results of the isotropic ladders 



^ - 1 (70) 



where we have used the data of table 7 together the value A2/ J = 0.504 p9| , pSj . 

In agreement with (^), Chakravarty has recently derived the exponential fall off of the gap 
with Hi using the 2+1 NLSM [|l5l. In his approach a spin ladder of width ni at zero temperature 
and periodic boundary conditions along the rungs, is equivalent to a Heisenberg plane of infinite 
extent at a finite temperature inversely proportional to n^. This allows the use of the 2+1 NLSM 
results to study ladder systems. In particular the expressions for the spin gap and correlation 
length are given by ( isotropic ladders J' = J), |jl5| 



nS ^Ps L \ I 1 



(71) 



/^(2+l) 



47rV2 5Zp<; ^-^P V Zc a) \^ ttV2S Zp^ 

where L and a are the width and lattice spacing of the ladder {L/a = ne). The exponential 
behaviour of eqs(^) and ( [71] ) agree in the classical limit S ^ oc provided we choose Zc and 



1 - 



Zp^ as the renormalization constants Zc{S,ni) and Zp^(S,n£) defined in 



This gives further 



support to the finite size correction propose in (^ 



In tables 7 and 8 we give the values of the gap and correlation length of the 4 and 6 legged 
ladders obtained using: i) numerical methods (Quantum Monte Carlo ||3^, ^ and DMRG @), 
ii) the NLSM in 2+1 ( |]r3|) and iii) the finite size correction of the NLSM results of [|r^. For 



the two later set of data the values of the renormalization constant for S=l/2 are choosen as 
Z.. = 0.724 |1[. 



1.1^ 





DMRG 


QMC 


NLSM{2 + 1) 


NLSM{2 + l)+finite size 


4 


0.190 


0.16 - 0.17 


0.268 


0.209 


6 


? 


0.055 - 0.05 


0.064 


0.050 



Table 7: Ladder's spin gap. 



lie 


DMRG 


QMC 


NLSM{2 + 1) 


NSLM{2 + l)+finite size 


4 


5-6 


10.3 


6.23 


7.37 


6 


? 


~ 30 


26.2 


31.6 



Table 8: Ladder's correlation length. 

We observe from tables 7 and 8 that the finite size modification of the Chakravarty formulas 
([n|) seems to give a rather good agreement with the numerical results. Further work needs to 
be done to settle this matter. 
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Limits of Applicability of the Ladder's Map 

We mentioned at the beginning of the construction of the map from the ladder into a unique 
NSLM field, that it would be valid for the intermediate coupling region {J' / J ~ 0{1)). We shall 
next explain this point in more detail. 

First of all let us consider the weak coupling region J' / J « 1. As shown in eqs. (|38D the 
masses of the higher modes of the ladder, at the linearized level, are of order V JJ'. On the 
other hand the mass generated non perturbatively is given by (|69D . A consistent truncation of 
the massive modes then requires that the mass of these modes should be larger than the mass 
generated non perturbatively, which leads to, 

j^^-Bn, « L (72) 

This equation implies that there is a lower critical value, {J'/J)c below which the truncation 
of the high energy modes, at least in the way it is done here, is not valid. In particular, in 
the weak coupling region the spin gap is approximately proportional to J' ( for = 2, A ~ 
0.41J')|^9|, |0|. This behaviour is not consistent with (|69|). On the other hand eq.(72) 



suggests that the range of applicability of our model is bigger as rii increases. 

Let us now consider the strong coupling regime {J'/ J » 1), and choose as an example the 
two legged ladder. For the columnar block, the map (^) reads. 



(73) 



Si{n)= ll(n) + S(-l)"$(n) 
which plugged into the ladder Hamiltonian leads, without making any approximation, to 

i^iaddcr = E Y + J (^l(^) + 1) - 252$(n) $(n + 1)) (74) 

If we now apply a gradient expansion in the fields 1 and we obtain the results given in 
(1^) (see also table 9), for the case ra^ = 2. It is interesting to compare (|7^ with the lattice 
NLSM Hamiltonian ([T^. We see that the gradient expansion and truncation of the \s field is 
crucial in order that (0) becomes a NLSM Hamiltonian. In the strong coupling regime, as can 



be seen from the mean field analysis of what is pertinent is to drop the term proportional 
to l(n)l(n + 1) in (|7^. If we do this approximation then (^) become the discrete NLSM 
Hamiltonian with the following identification of g and c, 

1 HF (75) 



fl'latt — 5' 

If we replace these parameters into the gap formula obtained in the strong coupling regime 
(0) we get 

A/J' = l- ^-^ + 0.296^^ (^j (76) 

If we replace in ( |7^ ) S*^ by S{S + 1) and we particularize to the case S' = 1/2, then we get 
the correct behaviour of the gap in perturbation theory to order J/J' |l30[| . 
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The previous discussion illustrates that in the strong coupling limit one has to be careful in 
making gradient expansions of the NLSM fields. A more detailed analysis is needed to clarify 
this matter. 



Haldane Phase in the 2-Legged Ladder 

An important question concerning even spin ladders is wether these systems are in a fundamen- 
tally new state or they are in a more familiar state, as for example the integer spin chains. This 
problem has been addressed by various authors arriving at different conclusions p5| , In this 
section we shall apply the NLSM techniques to clarify this issue, finding support for the view of 
4^ that the 2 legged ladder is in the same phase than the S = 1 chain ||4^. A more detailed 
discussion is delayed to the future. 

The rii = 2 ladder can be studied using columnar, rectangular and diagonal blocks, 
values of the corresponding NLSM parameters are given in table 9. 



The 



Block 


6 


c 




Columnar 





1 v\ 1/2 

2JS (l + i^) 




Rectangular 







V ^ 2jJ 


Diagonal 


AnS 


2J5(1 + ^) 


1 

s 



Table 9: NLSM parameters of the 2 legged ladder 

From table 9 we see that the parameters obtained using the diagonal blocks coincide with 
those of an effective spin chain with spin and exchange coupling given by. 



eff 



25, Jeff = 7(2J+ J' 



4 



(7^ 



This is not an accidental fact. An alternative way to arrive to (^) is to construct an effective 
model in terms of the spin 2S states formed out by symmetrization of two spin S states located 
in diagonal positions of the ladder ( see fig. 5) [Q, i.e. 



Seff(n) = Si(n) + S2(n-l) (79) 

The Hamiltonian governing this effective spins is given, to lowest order in perturbation theory 
by the chain Hamiltonian, 



eff 



j')ESeff(n) Seff(n+ 1) 



^0) 



The NLSM parameters corresponding to this model (^) are precisely the ones given in table 
9 for the diagonal blocks. Hence the 2 legged ladder is mapped into the same NLSM model 
as the S=l chain, which suggests that both systems are in the same phase. The relationship 
between the "standard" phase of the ladder, given by the RVB picture of ||3^, and that of the 
spin 1 chain appears, from the point of view of the NLSM, as the relation between two different 
types of blockings, namely the rectangular and the diagonal ones. This relation is given by. 



1 = 1 + $ = $ 
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where 1, $ ( resp. 1, $) are the NLSM fields associated to rectangular (resp. diagonal) blocks. 
Eg . (^I]) is a cannonical transformation, identical to (|29|) , which changes 9 from into AttS, 
leaving invariant the values of c and g. This poses a puzzle since the later parameters are 
different for rectangular and diagonal blocks, except in the case J' = 2 J. This discrepancy must 
be understood along the lines of our discussion about the different values of g for columnar and 
rectangular blocks, and very likely do not affect the conclusion about the equivalence between 
the 2 legged ladder and the spin 1 chain. 



Spin Ladders with Dimerization 

The main property of uniform spin ladders is that there are no phase transitions as one varies 
the ratio J'/ J. It is thus interesting to investigate the existence of new phases by enlarging 
the parameter space of the model. There are plenty of possibilites at hand: dimerization, 
frustration, spin deffects, etc. We shall review here the first one. There are also many possible 
types of dimerization in a ladder. We shall consider below dimerizations only along the legs. 
In this case we may distinguish between columnar and staggered dimerizations, for which the 
intraleg coupling constant is given by, 

, . _ J J(l + (-1)"7) (columnar) 
' \ J(l + (-l)"+"7) (staggered) ^ ^ 

The dimerization parameter 7 will be choosen to vary in the interval (1,-1) in order not to 
change the AF character of the legs. The rung coupling constant J' may be positive or negative. 



so all together there are 4 types of models. In references were studied the models with 
columnar dimerization and ferromagnetic rung coupling. If J' —00 the dimerized ladder 
becomes effectively a dimerized chain and one can apply the results known for chains. 

A chain with spin S and alternation parameter 7 can be mapped into a NLSM with 6 given 

by 1,0, 



9 = 27iS{l + -f) (83) 

The criteria that the NLSM models with 6 = 71 { mod 2tt) are masless ^ implies the 
existence of 2S critical points Indeed, as 7 varies from -1 to 1, the parameter 6, given in 



(P3|), passes 2S times through vr. If S=l/2 there is a single critical point corresponding to the 
non dimerized chain (i.e. 7 = 0). If S=l there should exists two critical points for 7c = ±1/2. 
Numerical computations show that there are indeed two critical points located at 7c ~ ±0.25 
Hence the NLSM predicts the existence of these critical points but it is not precise about 



their localization. The NLSM prediction has also been confirmed for S= 3/2 and S=2 ||50 |. 

Returning to ladders with columnar dimerization we expect, from the above discussion, that 
there should exist two critical lines in the plane (7, J' /J)- In the example of the spin 1/2 two 
legged ladder there should exist two of these lines emanating at J' / J = —00,7c = ±0.25 and 
ending at the origin (i.e. J' /J = 0,7 = 0) [^. This critical lines separate the Haldane phase 
associated to the strong ferromagnetic ladder and a dimerized phase of weakly coupled chains. 

The staggered dimerization with AF rung couplings has been studied in The behaviour 



of these ladders is very interesting and it is based on the map into the NLSM. The value of the 
9 parameter of a staggered ladder with AF rung coupling is given by [pl[] . 
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9 = 2nSn,{l + ^ U^) (84) 

This formula contains ( ^3l) as a particular case (for rii = 1). Using the properties of one 
can conjecture from (Q) the existence of 2Sni critical lines in the plane (7, J'/ J) In the 



particular case of the spin 1/2 two legged ladder there should exist two critical lines separating 
also the RVB phase of the uniform ladders and a dimerized phase of weakly couple chains (see 
figure 6). The pahse diagram of the 3 legged ladder (S=l/2) contains besides the uniform critical 
line (7 = 0) two more critical lines ending at the walls I7I = 1 (see figure 7). 

The two cases analysed above ( i.e. columnar /F-rung and staggered/AF-rung ) have a very 
similar phase diagram which suggests some kind of relationship between them . Other types 
of dimerizations will be considered in ||52||. 
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Figure Captions 



Fig.l.- A generic spin ladder with n£ legs. 

Fig. 2.- Blocking of the spin chain needed to define the NLSM variables out of the spin ones. 
Fig. 3.- Columnar Blocking (n^ = 4). 
Fig. 4.- Rectangular Blocking {rii = 3). 
Fig. 5.- Diagonal Blocking (n^ = 3). 

Fig. 6.- Phase diagrams of the 2 legged ladder with staggered dimerization for S=l/2 and S=l 
(inset). If I7I — 1 the ladders degenerate into the "snake" chains which are depicted in the 
margins of the figure. 

Fig. 7.- Phase diagram of the 3 legged ladder with columnar (lower half plane) and staggered 
dimerizations (upper half plane) . Observe the existence of 3 critical lines that emerge from the 
origin. 
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